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Abstract 

We discuss the atom-atom scattering problem across a Feshbach resonance in a 
two-dimensional dilute Bose gas at zero temperature, in the limit where the s- 
wave scattering length exceeds the width of the vertical confinement. We determine 
a tunable coupling-strength parameter and by controlling it we evaluate how the 
condensate wave function spreads out with increasing atom-atom repulsions. We 
also discuss the stability of the condensate in the magnetic-field regime where the 
coupling has become attractive. 
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1 Introduction 



An ultracold two-dimensional (2D) Bose gas exhibits in several respects a dif- 
ferent physical behaviour from its three-dimensional (3D) counterpart [1,2]. 
In the macroscopic limit phase fluctuations at any finite temperature destroy 
the long-range order associated with Bose-Einstein condensation, but still al- 
low a Kosterlitz-Thouless transition to a superffuid "quasi-condensate" state 
characterized by local phase coherence [3]. A true condensate is restored in 
a trapped 2D gas by modifications induced by the confinement in the den- 
sity of states [4] , but the magnitude and sign of the scattering amplitude for 
atom-atom collisions can depend on external system parameters [5]. 

In fact, when the scattering length becomes comparable with the trap size 
a new two-particle scattering problem needs to be solved. The T-matrix for 
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two-body collisions in vacuo at low momenta and energy, which is used in 
3D to obtain the collisional coupling parameter, vanishes in the strictly 2D 
limit as the axial width of the trapping potential becomes smaller than the 
s-wave scattering length [1,2]. It is then necessary to evaluate the scattering 
processes between pairs of Bose particles by taking into account the presence 
of the surrounding bosons through a many-body T-matrix formalism [6-8] . 



The 3D s-wave scattering length can be tuned by means of external electric 
or magnetic fields [9-11]. Within this process a resonance can appear at a 
particular value of the magnetic field and allow control of the coupling strength 
ranging from positive (repulsive) to negative (attractive) values. One may 
study new dynamical effects of collective oscillations in large condensates [12] 
and the dynamics of collapsing and exploding condensates [13]. Calculations 
of such Feshbach resonances have been given for 3D gases of Li [14], 23 Na [10], 
39 K and 41 K [15], 85 Rb and 87 Rb [11], and 133 Cs [9]. Experimental observations 
have been reported for 87 Rb [16], 85 Rb [13,17-19], 23 Na [20], and 133 Cs [21]. 
Feshbach resonances in quasi-2D atomic gases have been discussed by Wouters 
et al. [22], who have predicted a shift in the position of the resonance from 
the squeezing of the confinement. 



In this work we examine the tuning of two-body collisions by means of a 
magnetic field in a strictly 2D condensate inside a harmonic trap at zero tem- 
perature. By modifying the coupling strength we can control the dimensionless 
physical parameter rj which is the ratio between the mean-field interaction en- 
ergy and the harmonic-oscillator level spacing [23]. In the case of repulsive 
interactions (rj > 0), starting from a weakly interacting regime (rj <C 1) and 
up to a strongly interacting one (77 1), we evaluate the equilibrium den- 
sity profile as a function of rj by solving numerically a non-linear Schrodinger 
equation by a split-step Crank-Nicholson discretization scheme. In the case of 
attractive interactions (rj < 0) we can give an upper bound for the number 
of bosons which can condense without the condensate collapsing. At variance 
from the 3D gas (see for instance [24,25]), we have found that in 2D this 
critical number does not depend explicitly on the strength of the harmonic 
confinement. 



The paper is organized as follows. Section 2 summarizes the Feshbach reso- 
nance formalism and applies it to the 2D Bose gas. In Sec. 3 we evaluate the 
equilibrium density profiles for the case of repulsive interactions and discuss 
the conditions for a stable condensate in the attractive coupling regime. Some 
concluding remarks are presented in Sec. 4. 
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2 Essential formalism 



2.1 Overview on scattering theory 

We start by recalling some well-known facts about a collision between two 
particles with internal degrees of freedom governed by an internal Hamiltonian 

H int , 

H int \a) = e a \a) (1) 

where \a) and e a are the single-atom eigenstate and its corresponding eigen- 
value. For example the internal Hamiltonian can describe the hyperfine inter- 
action between nuclear and electronic spins and their Zeeman coupling to an 
external magnetic field B. A two-body collision is described by the Hamilto- 
nian 

2 2 

H =£^ + Y, H i nt + V(r) (2) 

in the center-of-mass system, with p the relative momentum and m r the re- 
duced mass (see for instance [26]). The Hamiltonian in Eq. (2) is the sum of 
a part with eigenstates \af3) ("channels") tending asymptotically to a sym- 
metrized product of separate-atom internal states \a) and and of a finite- 
range interaction V(r) which couples the channels. The energy associated with 
an eigenstate \a/3) is E a> p = e a + Sp + h 2 k^j 3 /2m r , where Kk a j3 is the relative 
momentum of the two incoming particles. The asymptotic scattering state 
^^(r) can be expanded on all scattering channels in the form 

<Mr) = exp(*k Q/3 ■ r) \a(3) + £ ^/(Kp, Kp f*^ ' ~ l«W , (3) 

a>,fl> r 

where the outgoing channel is denoted by \ct'f3') and the scattering amplitude 
faf ^ s proportional to the T-matrix element (k^/| T |k aj g). The asymptotic 

magnitude of the momentum for the \ol'(5') channel is k a ipi = \jlm r {E a p — e a r — 6/3')? 
from conservation of total momentum and energy in the scattering process. 
The channel is said to be open (closed) if the momentum takes positive real 
(imaginary) value. 

The whole Hilbert space describing the spatial and spin degrees of freedom is 
now divided into a subspace P for open channels and a subspace Q including 
all closed channels [27]. The state vector and the Schodinger equation 
H |\&) = E\^} are projected onto the two subspaces by means of projection 
operators P and Q satisfying the conditions P + Q = 1 and PQ = 0. The formal 
solution of the coupled equations obeyed by the two projected components is 

\* Q ) = (E-H QQ + i8)- 1 H QP \*p) (4) 
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and 

(E-H PP -H , PP )\*p) = 0, (5) 
where H PP = PHP etcetera, 

H' PP = H PQ (E - H QQ + iS^Hqp , (6) 

and 5 is a positive infinitesimal. The term H PP in Eq. (5) describes the Fes- 
hbach resonances. It represents an effective non-local, retarded interaction in 
the P subspace due to transitions to the Q subspace and back. The Hamilto- 
nian H PP + H' PP gives the effective atom- atom interaction in the open-channel 
subspace, the effects due to the existence of bound states being taken into ac- 
count through the term H PP . 



2.2 Feshbach resonances in the 2D coupling 



A Feshbach resonance results when true bound states belonging to the closed- 
channel subspace match the energy of open channels: transient transitions may 
then be possible during the collision process [10]. Setting H PP = H + Ui where 
H includes the kinetic energy of relative motion and the internal Hamiltonian, 
the scattering matrix elements in the P subspace between plane-wave states 
with relative momenta k and k' read 

(k'| T |k> ~ (k'| Tx |k> + (fi 1>+ k'| H' PP |ft li+ k) . (7) 

Here Ti = U± + U\GqTi represents the scattering process if Q subspace is 
neglected and = (1 — G £A) 1 is the wave operator due to the interaction 
potential U\, with Go = (E — Hq) -1 [28]. In Eq. (7) the effect of closed channels 
has been included at first order in H' PP . In the limit of zero relative velocity 
Eq. (7) leads to an expression for the coupling strength, 

where the sum is over all the states \ip n ) in Q subspace and we have set 
|f2i i+ k) ~ |ili i+ k') = \ipo). Here, g is the coupling strength estimated in the 
absence of closed channels and E t h is the threshold energy for the state |?/>o) 
at vanishing kinetic energy. Finally, when the threshold energy is close to the 
energy E res of one specific bound state \^ res ), this contribution will be domi- 
nant and the contributions from all other states may be included in g through 
an effective non- resonant scattering length a nr . For a strictly 2D condensate 
we get 

= 4irh 2 /m | (ipresl H QP \ifj ) | 2 
In \4h//j,mal r \ E th - E res 

where \i is the chemical potential of the condensate. 
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Equation (9) is our main result. The non- resonant first term has been cal- 
culated from the two-body T-matrix for a binary collision occurring at low 
momenta and energy in the presence of a condensate, which fixes the energy 
of each colliding atom at the chemical potential [29,30]. The resonant second 
term can be tuned by exploiting the dependence of the energy denominator on 
external parameters and in particular on the strength of the applied magnetic 
field. Explicitly, if the energy denominator vanishes at a value Bo of the field, 
then by expansion around this value we have 

E th - E res « (m res - m a - mp)(B - B ) (10) 

where m Q / g = —de a ^ / dB are the magnetic moments of the two colliding atoms 
in the open channel and m res = —dE res /dB is the magnetic moment of the 
molecular bound state. Equation (9) then yields 

- ( 1 AB 

m yln \AH/ iima 2 nr \ B — B 

where the width parameter AB is defined as 

AB = rn | (ip res \ H QP \ip } | 2 
Aith 2 m« s - m Q - tti/3 

Higher-order terms in H' pp will lead to a broadening of the resonance. However, 
the width of resonant states close to the threshold energy in the open channel 
is usually small because of the low density of states, and Feshbach resonances 
are usually very sharp [28]. 

In the next Section we will analyze the consequences of tuning the coupling 
strength on the equilibrium state of a 2D condensate in a harmonic trap. 



3 Equilibrium properties as functions of the coupling 

We consider a Bose gas with a fixed number N of particles in a 2D harmonic 
potential V ext (r) = mu 2 r 2 /2. At zero temperature the condensate wave func- 
tion \P(r) is described by the mean-field Gross-Pitaevskii equation 

tf(r) = fi* (r) , (13) 

where n(r) = ^(r)! 2 is the density of the condensate. We introduce the di- 
mensionless coordinate x = rja^ where = (h/muj) 1 / 2 is the harmonic 
oscillator length, and the dimensionless chemical potential // = 2/i/hcu. Thus 
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Eq. (13) can be written in rescaled form as 



d 2 



1 d 



dx 2 x dx 



V(x) +x 2 m(x) + 2r] 



*(0) 



= fi'V(x) , (14) 



where 77 = g\^(0)\ 2 /hu is a dimensionless coupling parameter given by the 
ratio between the mean-field interaction energy per particle and the oscillator 
level spacing. The normalization condition on the wave-function is 



roo 

N = a\ \ xdx\^(x)\ 2 . 
Jo 

Equations (14) and (15) need to be solved self-consistently. 



(15) 



3.1 Repulsive coupling (rj > 0) 



In the case of repulsive interactions Eq. (14) always admits a solution: the 
condensate exists and is stable in the 2D harmonic trap. Of course, far away 
from the resonance the dimensionless coupling parameter 77 is primarily deter- 
mined by the non-resonant scattering length. On approaching the resonance 
by increasing B towards B , the coupling parameter becomes very small for 
B ~ B — AB/(\n \ Ah/ fxma\ r \): the non-linear term in Eq. (14) becomes small 
and the condensate wave function is close to a Gaussian profile. On further 
increasing B the gas enters the strong-coupling regime and in the limit 77 ^> 1 
the kinetic energy term in Eq. (14) and the discrete structure of the trap levels 
become irrelevant. The profile then approaches a Thomas- Fermi form [31,32]. 
The shape of the condensate density profile, obtained by the numerical solu- 
tion of Eq. (14) at several values of the coupling parameter 77, is depicted in 
Fig. 1. The Figure illustrates comparisons with the free-gas and Thomas- Fermi 
profiles, according to the discussion given just above. 



3.2 Attractive coupling (rj < 0) 



On going through the resonance the Thomas-Fermi-like density profile will 
suddenly shrink and collapse. Indeed, in the strong-coupling attractive regime 
the kinetic energy does not suffice to compensate for the negative mean-field 
energy term and Eq. (14) does not admit physical solutions. A condensate 
with a limited number of particles can only exist in the weakly attractive 
regime. An upper bound for the critical number N c of bosons in a condensate 
with weakly attractive interactions can be estimated by comparing the kinetic 
energy and the mean-field energy at the centre of the trap. This leads to the 
inequality 

1 + 2?7 > (16) 
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Fig. 1. Scaled density profile |^ / (x)| 2 /|^ / (0)| 2 as a function of the dimensionless radial 
distance x = r ja\ lo for various values of the coupling parameter rj. Dashed curves 
are Thomas-Fermi profiles and the dotted curve is the free-gas profile. 

for condensate stability, where in the calculation of the kinetic energy we have 
used the free-gas wave function. Using the expression of rj we find from Eq. 
(16) 



1 AB 

+ 



(17) 

\n\4h/fima 2 nr \ B — B y J 

At variance from the trapped 3D gas [24], the critical number N c does not 
depend on the harmonic oscillator length, the dependence being cancelled by 
the fact that in 2D the mean-field energy scales as ~ % G 2 like the kinetic 
energy. Let us remark that the condensate in the case of attractive coupling 
can be metastable even if the number of condensed bosons is lower than N c . A 
three-body loss term should be included in Eq. (13) in order to calculate the 
decay time and to describe the time evolution of such a condensate, as was 
done for the 3D Bose gas [25,33,34] and very recently for a 3D boson-fermion 
mixture [35]. An evaluation of three-body losses in the 2D Bose gas and of the 
time evolution of a 2D condensate in the weakly attractive coupling regime 
may be considered in the future. 



4 Concluding remarks 



In summary, we have applied the formalism already developed to treat Fes- 
hbach resonances in 3D boson gases to the case of a 2D Bose gas at zero 
temperature. One of the main results is that an external magnetic field could 
be used to drive the coupling from repulsive to attractive even in a pancake 
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geometry as the scattering length becomes larger than the axial width of the 
condensed cloud. For the case of repulsive interactions we have evaluated the 
particle density profiles at varying coupling strength, demonstrating the ap- 
proach to a Thomas- Fermi profile followed by collapse. In the opposite regime 
of attractive coupling we have given an estimate for the critical number of 
bosons that can form a stable condensate. We have found that this critical 
number is inversely proportional to the absolute magnitude of the coupling, 
independently of the strength of the in-plane confinement. 
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